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ABSTRACT 

The starburst galaxy M82 shows a system of Ha-emitting filaments which extend to each side of 
the galactic disk. We model these filaments as the result of the interaction between the winds from a 
distribution of Super Stellar Clusters (SSCs). We first derive the condition necessary for producing a 
radiative interaction between the cluster winds (a condition which is met by the SSC distribution of 
M82). We then compute 3D simulations for SSC wind distributions which satisfy the condition for a 
radiative interaction, and also for distributions which do not satisfy this condition. We find that the 
highly radiative models, that result from the interaction of high metallicity cluster winds, produce a 
structure of Ha emitting filaments, which qualitatively agrees with the observations of the M82, while 
the non-radiative SSC wind interaction models do not produce filamentary structures. Therefore, our 
criterion for radiative interactions (which depends on the mass loss rate and the terminal velocity of 
the SSC winds, and the mean separation between SSCs) can be used to predict whether or not an 
observed galaxy should have associated Ha emitting filaments. 

Subject headings: Hydrodynamics - shock waves - star: winds, outflows - galaxies: star cluster - 
galaxies: starburst - galaxies: individual: (M82, NGC253) 



1. INTRODUCTION 

Supergalactic winds are a very important mechanism 
in the evolution of the Universe. Lynds & Sandage (1963) 
found a large outflow in M82 and gave the first defini- 
tion of a starburst. They proposed that material was 
ejected from the nuclear regions as the consequence of 
a SN explosion. Different analytical models have been 
put forward to explain the so-called supergalactic winds 
(SGWs), see, e. g., the work of Chevalier & Clegg (1985), 
Tomisaka & Ikeuchi (1988), Heckman et al (1990), Lei- 
therer & Heckman (1995), Strickland & Stevens (2000), 
Canto et al (2000), Tenorio-Tagle et al. (2003), and 
Cooper et al. (2008). In the nearest starbursts, they 
have shown that the supergalactic winds (SGW) form 
through the collective effects of many individual stellar 
cluster winds, which are in turn formed by the collective 
effect of many individual stellar winds. 

The best-studied starburst galaxy is M82. It has an 
extended, biconical filamentary structure in the optical 
(Shopbell & Bland-Hawthorn 1998; Ohyama et al. 2002). 
This optical emission is embedded in a pool of soft X-ray 
emission, detected with the Chandra X-Ray Observatory 
(Griffiths et al. 2000) as well as with XMM-Newton 
(Stevens et al. 2003). This X-ray emission extends sev- 
eral kpc away from the nuclear region, into the IGM. 

With the Hubble Space Telescope (HST), it has been 
possible to isolate the building blocks of some star- 
bursts:Super Stellar Clusters (SSCs). This objects are 
very massive and dense stellar clusters, with masses in 
excess of 10 4 M enclosed in radii of 3-10 pc. They are 
young, and typically contain from ~ 100 up to ~ 10 4 
OB stars. In M82, the SSCs Ha luminosities are in the 
range of (0.01-23) x 10 38 erg s" 1 (Melo et al 2005). These 
authors catalogued 197 SSCs in M82 with masses in the 
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10 4 < M/M < 10 6 range. This exceptional density of 

massive clusters (620 kpc 2 ) produces a suitable scenario 
for a study of the interaction between cluster winds. 

The differences in the populations of SSCs of several 
starburst galaxies can help us to understand which prop- 
erties are more important for producing filaments. For 
instance NGC253, which does not show a clear filmen- 
tary structure (Rice 1993), has a significantly lower SSC 
mass loss rate than M82. However, other authors claim 
that the lack filamentary structure in this galactic wind 
is simply a result of not having deep enough Ha observa- 
tions of this objects (Hoopers et al. 1995). In NGC253 
Melo (2005) catalogued a total of 48 SSCs with an av- 
erage distance between them of 31 pc, in contrast with 
the mean separation between SSCs of 12 pc found in 
M82. In M82, a rich structure of filaments, extending 
out from the galactic plane, is present. There are many 
other examples such as NGC1569 (Anders et al. 2004; 
Westmoquette, Smith & Gallagher III 2008) and M83 
(Harris et al. 2001), that have high density of SSCs with 
similar masses to the ones of M82. 

Tenorio-Tagle et al. (2003) showed that the interaction 
of stellar winds from a collection of energetic, neighbour- 
ing SSCs could form a filamentary structure similar to 
the one observed in M82. These authors propose the for- 
mation of such filaments by means of supernovae explo- 
sions that would expel a huge amount of heavy elements 
into the Interstellar Medium (ISM), thus enhancing the 
radiative cooling of the outflows. 

In their model, Tenorio-Tagle et al. (2003) explained 
the SGW structure as a consequence of the interac- 
tion of winds from very close SSCs, in which stationary 
oblique shocks are responsible for shaping the material 
into dense and cold, kiloparsec-sized filaments. More re- 
cently, Cooper et al. (2008) presented a numerical study 
of a starburst-driven galactic wind. Their setup consisted 
of a series of neutral dense clouds placed in the galaxy 
disk. The dense clouds are swept up by the main shock 
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TABLE 1 

Cooling distance as a function of shock velocity 



v s [km s 1 ] 



logio(d c ,i [pc]) 



100 


-1.468 


-2.010 


-2.283 


200 


-1.194 


-1.960 


-2.259 


300 


-0.657 


-1.439 


-1.796 


400 


0.595 


-0.091 


-0.399 


500 


0.957 


0.269 


-0.037 


600 


1.176 


0.516 


0.209 


700 


1.449 


0.763 


0.455 


800 


1.601 


0.913 


0.605 


900 


1.727 


1.038 


0.714 


1000 


1.893 


1.208 


0.900 


1100 


2.108 


1.438 


1.123 



wave, leaving behind a filamentary structure. 

In the present paper, we discuss models for the for- 
mation of filamentary structures as the result of the in- 
teraction between winds from a system of many clusters 
with a disk-like spatial distribution. In these models we 
consider only the flow resulting from the wind interac- 
tions, and do not include the effect of the stratified ISM 
present in the region of the galactic disk. 

The paper is organized as follows. In §2, we study 
which combinations of parameters (mass loss rate M w , 
cluster wind velocity v w , and separation between nearby 
star clusters D) yield a highly radiative SGW flow. In 
§3 we describe a series of 3D simulations of the inter- 
action of cluster winds in this highly radiative regime. 
We present predictions of the emission in the optical and 
X-ray wavelength ranges, and compare them with obser- 
vations of galactic winds (i.e. M82, NGC253, etc.) in §4. 
Finally, in §5 we present our conclusions. 

2. THE FORMATION OF FILAMENTARY STRUCTURES 
2.1. General considerations 

The collective effect of the individual stellar winds in- 
side an SSC result in a cluster wind. Outside the outer 
radius of the cluster (i. e., for radii r > r c , where r c 
is the cluster radius), the cluster wind has an approxi- 
mately constant velocity, and a density that falls oc r~ 2 
(where r is the distance to the cluster centre, see, e. g., 
Rodriguez- Gonzalez et al. 2007). 

Cold filaments in SGWs can be formed by the the in- 
teraction of such SSC winds, provided that the cooling 
is efficient. Since the gas density of the cluster wind falls 
with distance, this is likely to happen only when the clus- 
ters are very close each other. The terminal velocity of 
the SSC winds and their mass loss rate are also important 
for determining whether or not we have efficient cooling. 

2.2. Radiative losses in a cluster wind 

Let us consider a galaxy with a central stellar cluster 
density n (number of clusters per unit volume). For the 
sake of simplicity, we consider stellar clusters with iden- 
tical, isotropic winds, mass loss rate M w , and a terminal 
wind velocity v w . The typical separation between stellar 
clusters is D w n -1 / 3 . 

Therefore, two- wind shock interactions between nearby 
clusters pair will occur at a typical distance ~ D/2 from 
each of the stellar clusters, so that the pre-shock densities 
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Fig. 1. — This diagram shows M w x (10 pc/D) as a function of 
v w (where D is the separation between clusters, M w is the mass 
loss rate, and v w the cluster wind velocity) for a "cooling param- 
eter" = d coo i/D = 1. Three curves are shown, corresponding to 
models with 1, 5, and 10 times solar metallicities. The region above 
the curves represents the parameter space in which the two-wind 
interactions of cluster winds are radiative. 



would have values of 



v pre 



1.2>mH^D 2 v u 



(1) 



where tuh is the Hydrogen mass, and we have assumed 
a 90% H and 10% He particle abundance. 

The shock interactions between nearby clusters will be 
radiative if the cooling distance d coo \ satisfies the condi- 
tion 

d cool ^ 1 ^ 



K = 



D 



< 1. 



In order to estimate the cooling distance, we proceed 
as follows. Consider the structure that will be formed 
by the leading shock of a star cluster wind with velocity 
« v w (i. e. the cluster wind velocity) with a preshock 
density given by equation (1). For a cooling function in 
the low density regime, the cooling distance scales as the 
inverse of the pre-shock density, yielding: 



d>cool{P j prei r V s) — 



d Ci i(v s ) 



°pre 



(3) 



where v s (= v w ) is the shock velocity, and d Ci i(v s ) is the 
cooling distance (to 10 4 K) behind a shock with n pre = 
1 cm -3 . We have calculated the function d Cj i(v s ) by 
integrating the equations for a ID, stationary shock (see, 
e. g., Shull & McKee 1979) using the coronal equilibrium 
cooling functions for 1, 5 and 10 times solar metallicities 
of Raymond, Cox & Smith (1976). The resulting values 
of the cooling distance d Ci i(v s ) (to 10 4 K, for a n pre = 
1 cm -3 preshock density) are given in Table 1. 

Setting k = 1 in equation (2) and also using equation 
(3), one obtains 



M 
~D 



= 1.3m H 7rv w d cA (v w ) . 



(4) 



which gives the relation between M, D and v w for SSCs 
that are just entering the highly radiative regime. In 
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Figure 1, we plot the ^ ratio as a function of v w and for 
SSCs with k = 1 (see equation 4) for the three metallic- 
ities which we have considered (see above and Table 1). 
The curves that would be obtained for n < 1 always lie 
above the n = 1 curve corresponding to the same metal- 
licity. Figure 1 shows that in a substantial part of the pa- 
rameter space which would be expected for galaxies with 
low or intermediate mass young stellar clusters, cooling 
parameters n < 1 are produced. A galactic wind in this 
"highly radiative regime" will have dense, cool structures 
resulting from the radiative shocks in the interactions be- 
tween nearby stellar clusters. In the following section, we 
present numerical simulations of flows in this regime. 

3. NUMERICAL SIMULATION 

We have computed 3D numerical simulations with the 
full, radiative gasdynamic equations. The simulations 
solve a multiple stellar wind interaction problem with the 
3D, adaptive grid "yguazu-a" code, which is described in 
detail by Raga et al. (2000, 2002). The first four sim- 
ulations were performed on a five-level binary adaptive 
grid with a maximum resolution of 0.488 pc (correspond- 
ing to 512x512x1024 grid points at the maximum grid 
resolution), in a computational domain of 250x250x500 
pc. We have also computed a fifth simulation also on a 
five-level adaptive grid, but with a maximum resolution 
of 0.976 pc, corresponding to 1024 x 1024 x 2048 grid 
points (at the maximum grid resolution) in a domain of 
1 x 1 x 2 kpc (see Table 2). 

In order to study the formation of filaments in the ra- 
diative regime we ran 6 models. We have assumed that 
the computational domain was filled initially with a ho- 
mogeneous, stationary ambient medium with tempera- 
ture T env = 1000 K and density n env = 0.1 cm -3 . This 
environment clearly does not include the ISM of the disk 
of a starburst galaxy, for which densities of ~ 1000 cm -3 
might be more appropriate. We choose this low density 
so that the ISM in the region in between the cluster is 
rapidly flushed away by the cluster winds, and a station- 
ary wind interaction structure is reached as rapidly as 
possible. If we introduced a higher ISM density (within 
the galactic plane), there would be a transient regime 
(lasting longer for higher ISM densities) before the clus- 
ter winds break out of the galactic plane and travel into 
the intergalactic medium. The density which we have 
chosen is closer to a possible density for the inter galactic 
medium into which the cluster winds will be expanding 
once they leave the galaxy. 

We model the SSCs as spherical wind sources, centered 
at the positions as described below. For models M1-M4 
(see Table 2), the wind sources have a radius of R c = 2.92 
pc (corresponding to 6 pixels at the maximum resolution 
of the adaptive grid, which is always present at the wind 
sources), for models M2b and M3b their radii is R c = 
5.84 pc (also corresponding to 6 pixels at the highest 
resolution grid). 

Within these spheres, we impose (at every timestep) a 
wind with a temperature of T c = 10 7 K, an outwardly 
directed velocity v c = 1000 km s -1 , and a density that 
follows an r~ 2 law (where r is the radial coordinate mea- 
sured outwards from the position of each wind source), 
scaled so as to obtain the correct the mass loss rate. 

The positions of the wind sources are obtained by sam- 
pling a uniform random distribution. However, in prac- 



TABLE 2 
Model parameters 



Model 


Number of 


z 


Dt 


M S c 




SSCs 


[Zq] 


[pc] 


[M© yr- 


Ml 


100 


1 


9.27 


2 x 10- 


Mia 


100 


1 


9.27 


2 x 10- 


M2* 


100 


5 


9.27 


2 x 10- 


M3* 


100 


10 


9.27 


2 x 10- 


M4 


15 


10 


21.88 


4 x 10- 


M4a 


15 


10 


21.88 


4 x 10- 



TThe value quoted here is the mean distance between neighbouring 
clusters (nearest neighbors) as measured from the positions that result 
from sampling a uniform random distribution. 

■^We have also run models M2b and M3b with the same parameters 
as M2 and M3, but with 1/2 of the spatial resolution and 4 times the 
spatial extent as M2 and M3. 

tice, we have to modify the obtained cluster distribu- 
tions to avoid the overlap of sources. The stellar cluster 
sources are placed in a cylindrical structure, meant to 
model the central region of a galactic disk. 

3.1. The starburst regions 

In order to reproduce the properties of M82 and of 
NGC253, we use the average properties reported by Melo 
et al. (2005) and Melo (2005). For M82, the mean 
stellar cluster mass (Msc) is of 1.7 x 10 5 M©, and the 
mean separation between neighbouring clusters_(A) is of 
12 pc. For NGC253, M 5C =6.2xl0 4 M and A=31 pc. 
We should note that the clusters of Melo et al. (2005) 
were detected at optical wavelengths. As there is a high 
level of obscuration in the nuclear regions of M82 as well 
as NGC253, the real number of clusters will be larger, 
and the mean separation between clusters will be smaller 
than the values quoted above. 

We estimate the mass loss rate Msc using the star- 
burst 99 models (Heckman & Leitherer 1995, Leitherer 
et al. 1999 2 ) with the appropriate cluster masses and 
a Salpeter initial mass function (IMF) including stars 
between 1 and 100 M©. The resulting mass loss rate 
(through winds and supernovae) is reasonably constant 
after ~ 4 Myr of the onset of the starburst, with some en- 
hancement when input by SNe overcomes the input from 
winds. This enhancement occurs at around t ~ 10 Myr, 
and lasts for about 10 Myr, then the mass loss rate re- 
turns to a value similar to the pre-supernovae phase for 
another ~ 20 Myr. Finally, the mass loss rate drops dra- 
matically once the last star of ~ 8 M© explodes as a 
supernova leaving only intermediate and low-mass star 
winds (at this point the star cluster would no longer be 
classified as SSC). 

For models M1-M4 we used a constant mass loss rate 
with a value consistent with the average mechanical lu- 
minosity of the SSCs. However, given the assumptions 
we have made (of having a system of identical clusters 
with average properties) a larger mass loss rate could be 
possible, especially if the SSCs are in the SN phase. Such 
a large mass loss rate could have a significant effect in the 
production of filaments. For this reason we have run two 
additional models, Mia and M4a, with extreme values of 
Msc, which correspond to upper limits (within the un- 
certainties), at the peak of the mass loss rate durinng the 

2 http:/ /www. stsci.edu/science/starburst99/ 
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Fig. 2. — Column density maps of models M1-M4. The bar at the bottom right gives the logarithmic gray-scale (color-scale in the online 
version) in cm -2 . The images were obtained by integrating density along lines of sight parallel to the x-axis. The distributions correspond 
to time-integrations of (2.5,2.5,2.5,5) x 10 5 yr (for models M1-M4, respectively). 



SN phase. Aside of the increased mass loss rate, these 
two models are identical to Ml and M4, respectively (see 
Table 2). 

3.2. The models 

We computed four numerical simulations of SGWs 
driven by stellar clusters, with different metallicities, 
mass input rates, and number of clusters. All of the 
models were computed with identical stellar clusters with 
a cluster wind velocity v w = 1000 km s _1 , placed in a 



cylinder with a radius of 100 pc and a height of 20 pc. 
This cylinder represents the central part of the starburst 
in a galaxy as M82, with dimensions approaching those 
of an individual starburst clump (O'Connell et al. 1995, 
Westmoquette et al. 2007). In reality the starburst 
of M82 has spatial extent of ~ 500 pc (see O'Connel 
& Mangano 1978). Numerical simulations with the full 
spatial extent for the cluster distribution would be much 
more demanding computationally. From SB99, using the 
mean cluster mass for M82 (Melo 2005), we estimate a 
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Fig. 3. — Column density maps for the models with mass loss rate in the SN phase (Mia and M4a). The bar at the right gives the 
logarithmic gray-scale (color-scale in the online version) in cm -2 . The images were obtained by integrating the density along lines of sight 
parallel to the x-axis. The corresponding integration times are of (2.5,5) x 10 5 yr (for models Mia and M4a, respectively). 



mechanical luminosity of 6 x 10 39 erg s 1 per SSC. Us- 
ing Esc = Msc v w/2 relation, this luminosity translates 
into a cluster mass loss rate of ~ 2 x 10 -2 M yr -1 . This 
cluster mass loss rate has been used for computing mod- 
els Ml, M2, and M3. These three models have different 
metallicities : Z = 1, 5 and 1OZ (for Ml, M2 and M3 
respectively) . 

These metallicities have been chosen as representative 
values of the metallicities of winds from starbursts dur- 
ing the first 20 Myr time-evolution. Using the metal 
yields of Meynet & Maeder (2002), Tenorio-Tagle et al. 
(2005) have shown that the metallicity of the combined 
stellar winds remains at a value of ~ 0.5 Zq up to a 
t = 3 Myr evolutionary time, rapidly grows to ~ 15 Z Q 
at t ~ 7 Myr, and then gradually decreases to ~ 3 Z@ at 
t w 20 Myr. 

Similarly, for NGC253 we obtain a mechanical luminos- 
ity of 1.9 x 10 39 erg s _1 , which corresponds to a cluster 
mass loss rate ~ 4 x 10 -3 M yr -1 . This cluster mass 
loss rate was used for computing model M4. 

For models M1-M3, we have placed 100 clusters within 
the cylindrical volume described in the beginning of this 
subsection, which results in a mean separation between 
clusters l) = 9.3 pc, which is similar to the mean sepa- 
ration A = 12 pc between the clusters of M82 (see §3.1). 
Model M4 has a Z = 10Z@ metallicity, and has 15 clus- 
ters within the cylindrical volume, resulting in a mean 
separation between clusters of D = 21.9 pc, which is 
similar to the mean separation A = 31 pc between the 
clusters of NGC253 (see §3.1). The parameters of models 
M1-M4 are summarized in Table 2. 

With these parameters, models M2 and M3 lie above 
the n = 1 curves (for the appropriate metallicities), so 



that they are clearly in the highly radiative regime (see 
Figure 1). Model Ml and M4 lie below the n = 1 curves, 
so that they are not in this regime. Therefore, we would 
expect to form dense filaments only in models M2 and 
M3. We evolve models Ml, M2, and M3 up to an integra- 
tion time of t = 2.5 x 10 5 yr, and M4 up to t = 5 x 10 5 yr. 
In Figure 2, we show the column densities obtained by 
integrating the number density along the x-axis of the 
domain, for the flow stratifications obtained at the end 
of the time-integrations. From this Figure, it is clear that 
models M2 and M3 produce a structure of ~ 10 — > 20 
filaments both above and below the plane of the galactic 
disk (within which are distributed the SSCs). Models 
Ml and M4 do not produce filamentary structures. 

The bipolar distribution of the filaments is a direct re- 
sult of the flattened cylindrical distribution of the SSCs 
(see Figure 2). We have also done simulations in which 
the cylindrical distribution has a height equal to the di- 
ameter (not included in the present paper), and they 
show filamentary structures with an almost isotropic dis- 
tribution. The same result is of course also found for a 
spherically symmetric cluster distribution. 

In Figure 3, we present the column density obtained 
in the models with extreme mass loss rates (Mia and 
M4a, see Table 2) at t=2.5 and 5 xlO 5 yr, respectively. 
Due to the increase of M, the cooling distances become 
an order of magnitude shorter (see the density scaling in 
equation 3), thus the adimensional parameter n for Mia 
and M4a is comparable to that of M3 (ft < 1). Therefore, 
as expected, Figure 3 shows filamentary scructure, even 
when it is absent for the same configuration at a lower 
mass loss. 
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Fig. 4. — Ha maps obtained for models M1-M4. The bar at the bottom right, gives the logarithmic gray-scale (color-scale in the online 
version) in erg s — 1 cm -2 sterad - 1 . The images were obtained by integrating the emission coefficient along lines of sight parallel to the 
x-axis. The distributions correspond to time-integrations of (2.5,2.5,2.5,5) x 10 5 yr (for models M1-M4, respectively). 

has lower radiative losses due to its lower metallicity), 
the Ha filaments appear only at a distance of ~ ±50 pc 
from the galactic plane. Models Ml and M4 do not pro- 
duce Ha-emitting filaments. The Ha filaments in M82 
(Ohyama et al. 2002) appear to extend down to the 
galactic plane, though this effect could partially be due 
to the fact that the plane of M82 is at an angle of ~ 10° 
with respect to the line of sight. We would therefore con- 
clude that the Z = 10Z Q metallicity model M3 appears 
to be more appropriate for reproducing the Ha filaments 



4. THE GALACTIC WIND EMISSION 
4.1. The Ha emission 

In Figure 4, we show the Ha maps computed by in- 
tegrating the Ha emission coefficient along the x-axis. 
The emission coefficient is computed with the interpola- 
tion formula given by Aller (1987) for the temperature 
dependence of the recombination cascade. 

In model M3, we see a structure of Ha filaments which 
extend out from the galactic plane. In model M2 (which 
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pc 

Fig. 5. — X-ray maps for models Ml-M4.The bar at the bottom right, gives the logarithmic gray-scale (color-scale in the online version) 
in erg s — 1 cm -2 sterad - 1 . The images were obtained by integrating the X-ray emission coefficient along lines of sight (parallel to x-axis), 
in the energy band of 0.3 to 2 keV. The distributions correspond to time-integrations of (2.5,2.5,2.5,5) x 10 5 yr (for models M1-M4, 
respectively). 

of M82. 



4.2. The X-Ray emission 

We have taken the density and temperature stratifica- 
tions from the SGW flow configurations, and used them 
to calculate the X-ray emission. We have done this by 
computing the emission coefficient in the 0.3 — > 2 keV 



photon energy range using the CHIANTI 3 atomic data 
base and software (see Dere et al. 1997, Landi et al. 
2006). For this calculation, it is assumed that the ion- 
ization state of the gas corresponds to coronal ionization 

3 The CHIANTI database and associated IDL pro- 
cedures, now distributed as version 5.1, are freely 
available at: http://wwwsolar.nrl.navy.mil/chianti.html, 
http:/ /www. arcetri.astro.it/science/chianti/chianti. html, and 
http: / / www.damtp.cam.ac.uk/user/astro/chianti/chianti.html 
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(a) M2b t=K2x10 6 yr (b) M3b t=1.6x10 6 yr 



500 



-500 




5x10" 



5x10" 



5x10" 



"O 

o 



E 

o 

i 

0) 

CP 
l_ 

(D 



7x10" 



5x10" 6 



-400 -200 
pc 



200 400 



Fig. 6. — Ha maps obtained for models M2b (left) and M3b (right). The bar on the right, gives the logarithmic gray-scale (color-scale in 
the online version) in erg s _1 cm -2 sterad -1 . The images were obtained by integrating the emission coefficient along lines of sight parallel 
to the x-axis. The distributions correspond to time-integrations of (1.2, 1.6) x 10 6 yr (for models M2b and M3b, respectively). 



equilibrium in the low density regime (i. e. the emission 
coefficient is proportional to the square of the density). 

The X-ray maps corresponding to the end of the time- 
integrations of models M1-M4 are shown in Figure 5. 
These maps show that the models Ml and M4 (with 
low radiative losses) produce an extended, diffuse X-ray 
emission. On the other hand, the highly radiative models 
M2 and M3 produce X-ray emission concentrated to the 
galactic plane. 

M82 shows X-ray emission which extends to large dis- 
tances (few kpc) from the galactic plane (Stevens et al. 
2003) . Because the preferred model for this galaxy is M3 
(which produces Ha filaments which start from the galac- 
tic plane, see Figure 4 and §4.1), we see that in order to 
reconcile the X-ray map predicted from this model (see 
Figure 4) and the observed X-ray distribution of M82 we 
need to assume that the observed X-ray emission comes 
from another component. This component, e. g., could 
be the shock of the expanding SGW against the inter- 
galactic medium. The emission from this shock is seen 
in the top and bottom parts of the frames corresponding 
to models M2 and M3 in Figure 5, but the shock struc- 
ture is already escaping from the computational domain. 

4.3. The Ha emission at larger distances from the 
galaxy 

Models M1-M4 extend only to 250 pc on either side 
of the galactic plane. As the filaments in M82 and 
NGC1569 extend to ~ 0.5 — > 1 kpc away from the plane 
of the galaxy, we have run models M2b and M3b (with 
the same parameters as M2 and M3, respectively, but 
with lower resolution and larger spatial extent, see §3 
and Table 2) in order to explore whether or not the Ha 



filaments predicted from our models do extend out to 
~ 1 kpc. 

In Figure 6, we show the Ha maps predicted from mod- 
els M2b (at t = 1.2 xlO 6 yr) and M3b (at t = 1.6xl0 6 yr). 
The two times are chosen so as to show an approximate 
time sequence of the flow, as the two models produce 
similar results. It is clear that in the two models, the 
Ha filaments extend to distances of ~ ±500 pc from 
the galactic plane. In the Ha map of M3b, there is a 
transition to a more complex morphology at distances 
> ±300 pc from the galactic plane, which is not seen in 
model M2b (in which basically the same filaments are 
seen at larger distances). 

5. CONCLUSIONS 

We study the formation of filaments in the galactic 
winds driven by young SSCs in regions of high star for- 
mation. We study models in which the filaments are pro- 
duced solely by the interaction between winds from SSCs, 
and do not consider the possible effects of the stratified 
or clumped ISM present in the galactic disk. In our mod- 
els, the winds propagate into a low density, homogeneous 
environment, which could correspond to the intergalactic 
medium. 

We have shown how a (M w /D) vs. v w diagram can be 
used to determine if the interaction between the winds 
of SSCs are radiative, in which case cold filaments can 
be formed. With the (M w /D) vs. v w diagram (Figure 1) 
one can diagnose if a galaxy with a high star formation 
rate will show filaments in optical emission lines, or only 
a galactic wind with X-ray emission. At the same time, 
given an observation of a galactic wind with a rich fila- 
mentary structure, one could predict the number of SSCs 
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which are contained within the observed galaxy. 

We computed three models of dense starbursts (MI- 
MS) based on the observed parameters of M82, and a 
fourth model (M4), based on NGC253. The models are 
based on cluster distributions with mean separations be- 
tween clusters similar to the observed values, but with 
spatial extents which are considerably smaller than the 
total extent of the starburst regions. The distributions 
used in the model could correspond to individual "star- 
burst clumps" such as the ones observed in M82 by West- 
moquette et al. (2007). 

For the parameters of M82 we computed four models 
(M1-M3, and Mia) with different metallicities (from so- 
lar to 10 times solar). For the parameters of NGC253 
we computed two models (M4 and M4a) with a 10 times 
solar metallicity. Models M1-M4 asssumed a mass loss 
rate consistent with the average mechanical luminosity 
of the SSCs. Models Mia and M4a are obtained with 
rather extreme mass loss rates, only achievable during a 
short-lived phase in which the input from SN dominates 
over stellar winds. 

For the models based on M82, we see the formation 
of cold, Ha-emitting, filaments that extend out from the 
plane of the the galactic disk in the models with metallic- 
ities > 5 Z Q (models M2 and M3). We also see formation 
of dense filaments for the models with extreme mass loss 
rate at the supernova phase of the clusters (models Mai 
and M4a), where a more conservative mass loss rate did 
not yielded filaments (Ml and M4). 

For the model based on NGC253 (M4), no filaments 
are produced, even for a 10 Zq metallicity. This is con- 
sistent with the present observations of this object, in 
which no Ha-emitting filaments have been detected. The 
prediction obtained from our models would be that with 
the distribution of SSCs of NGC253 (in pre or post SN 
phase), no Ha filaments should exist. However, we have 
shown that during the SN phase (model M4a and sec- 
tion 3.1) the formation of cold filaments is possible in 
such high metallic winds. Hoopes et al. (2005) claim 
that the lack of evidence of filamentary structure in the 
galactic wind of NGC253 could simply be the result of 
not having deep enough Ha maps of this object. In the 
future we will see whether or not this result is confirmed 
by deeper observations. 

We have also computed X-ray maps in the energy range 
of 0.3 to 2 keV. For the highly radiative models (M2 and 



M3), we find that the X-ray emission is concentrated in 
a region close to the galactic plane. In order to reconcile 
this result with the observation of extended X-ray emis- 
sion in M82, we would have to assume that the emission 
comes from another component in the flow, which could 
be the shock of SGW against the intergalactic medium. 

Finally, we have run two models (M2b and M3b) with 
the parameters based on M82 (and abundances of 5 Z Q 
and 10 Z , respectively) with lower computational reso- 
lutions but with larger spatial extents. From these runs, 
we find that the models produce a filamentary Ha mor- 
phology that extends out to ~ 0.5 — > 1 kpc on each side 
of the galactic plane. 

We end by noting that in this paper we prove that the 
radiative interaction between the winds from a disk-like 
distribution of identical SSCs (super stellar clusters) does 
lead to the production of dense filaments extending out 
to 100 — > 1000 pc away from the galactic plane. These 
structures produce Ha emission that might correspond 
to the filaments observed in starburst galaxies such as 
M82 and NGC1569. Of course, in our models we do not 
consider many elements that are likely to be important 
in starburst galaxies. Important among these might be: 

• the presence of a dense, stratified galactic disk in- 
terstellar medium, 

• dumpiness in this medium, 

• the existence of a distribution of sizes and ages for 
the SSCs, 

• the effect of galactic mass loading. 

There is therefore a large field for future theoretical stud- 
ies on the formation of filamentary strutures in starburst 
galaxies. 
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